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Abstract 

An alternative derivation of the first-order relativistic contribution to perihelic precession is pre- 
sented. Orbital motion in the Schwarzschild geometry is considered in the Keplerian limit, and the 
orbit equation is derived for approximately elliptical motion. The method of solution makes use of 
coordinate transformations and the correspondence principle, rather than the standard perturba- 
tive approach. The form of the resulting orbit equation is similar to that derived from Newtonian 
mechanics and includes first-order corrections to Kepler's orbits due to general relativity. The as- 
sociated relativistic contribution to perihelic precession agrees with established first-order results. 
The reduced radius for the circular orbit is in agreement to first-order with that calculated from 
the Schwarzschild effective potential. The method of solution is understandable by undergraduate 
students. 
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I. INTRODUCTION 



The derivation of the relativistic contribution to perihelic precession is found in almost 
every general relativity textbook [1, 2, 3, 4, 5, 6, 7, 8, 9, 10]. The standard approach 
is perturbative and results in a first-order prediction for the relativistic contribution to 
perihelic precession which is in agreement with observations. In particular, the agreement 
of this contribution with the observed shift of the perihelion for Mercury is one of the 
most celebrated verifications of Einstein's general theory of relativity [11, 12, 13, 14], and 
it is important that students are presented with a convincing calculation of this result. We 
present an alternative approach to this problem, and obtain an approximate relativistic orbit 
equation. Perihelic precession arises as one of several small corrections to Kepler's orbits, 
and agrees with established first-order calculations. 



The path of a small test mass near a spherically symmetric central mass M is uniquely 
described by the Schwarzschild line element [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 15, 16, 17, 18, 19]: 



where d(2 2 = d9 2 + sin 2 9 dip 2 is an infinitesimal element of solid angle, and r s = 2GM/c 2 
is the Schwarzschild radius. (G is Newton's universal gravitational constant, and c is the 
speed of light in vacuum.) We parameterize timelike geodesies as ds 2 = c 2 dr 2 , where r is 
the proper time along the path of a test particle. The variational principle for a geodesic 
results in Lagrange's equations, 



ds 2 = (1 - r s /r)c 2 dt 2 - (1 - r s /r) -1 dr 2 - r 2 dJ? 2 , 



(1) 
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dr dx^ 
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for each of {x^} = {t, r, 9, </?}, where x^ 



dx^/dr, and 



£ = (1 - r s /r)c 2 i 2 - (1 - rJryH 2 - r 2 {9 2 + sin 2 9<p 2 ). 



(3) 



Consider orbits in the plane defined by 9 = tt/2, so that 9 = 0. The invariants of the 
motion are given by Lagrange's equations for ip and t: 




which implies that £ = r 2 (p = constant; and 



2c 2 i(l-r s /r)i] = 0, 



(5) 
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which implies that k = (1 — r g /r) i = constant. The first of these invariants is the relativistic 
analogue to the Newtonian equation for the conservation of angular momentum per unit 
mass. The radial equation is most easily obtained by dividing the line element Eq. (1) by 
dr 2 , resulting in 

c 2 = (1 - r s /r)c 2 i 2 - (1 - rg/r^f 2 - r 2 ip 2 . (6) 

In terms of the invariants of motion, and with the definition E = \{k 2 — l)c 2 , the radial 
equation Eq. (6) may be expressed as 

^ i. 2 gm e £% 

E =2 r - — + 2^-2^- (7) 



Time is eliminated by the replacement 

nr — — / 

dip r ' 



r = -£^--, (8) 
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resulting in 

F! = 

2 \dip r J r 2r 2 2r 3 
Differentiating Eq. (9) with respect to ip results in a nonlinear differential equation for 

= f — --GM + -- 3 ^ (10) 

We anticipate a solution of Eq. (10) that is near Keplerian and introduce the radius of a 
circular orbit for a classical particle with the same angular momentum, r c = £ 2 /GM. The 
result is 



d 2 r c r r . 



+ -± = 1 + 3e 



d(y9 2 r r 

where e = (GM/ic) 2 = r s /2r c . The conic-sections of Newtonian mechanics [20, 21, 22] are 
recovered by setting e = 0: 

' |2 ''" + = 1, (12) 



dip 2 r r 
which implies that 

T 

— — 1 + ecosy?, (13) 
r 

where e is the eccentricity. 
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II. KEPLERIAN LIMIT AND PERIHELIC PRECESSION 



Relativistic orbits near a spherically symmetric massive object, such as our Sun, are 
described by Eq. (11). The planets of our solar system are described by near circular orbits. 
Mercury has the largest eccentricity (e pa 0.2), and the next largest is that of Mars (e pa 0.09). 
If e is taken to be a small relativistic correction to the near-circular orbits of Newtonian 
mechanics, it is convenient to make the change of variable 1/er = r c /r — 1 <gC 1. The last 
term on the right-hand-side of Eq. (11) may now be approximated as (r c /r) 2 pa 1 + 2/cr, 
resulting in a linear differential equation for l/a(ip): 

1 d 2 1 1 - 6e 1 



3e dif z a 6e a 
The additional change of variable a = ipy/l — 6e results in the familiar form: 

d2 °"c , 0"c - , 1K v 

dcr a a 

where a c = (1 — 6e)/3e. The solution is similar to that of Eq. (12): 

— pa 1 + A cos a, (16) 

(7 

where A is an arbitrary constant of integration. In terms of the original coordinates, Eq. (16) 
becomes 

T 

— pa 1 + ecos kip, (17) 
r 

where 

1 - 6e 



L l-3e' 
3eA 



(18) 

-l-fc' (19) 
K=(l-6e)3. (20) 

According to the correspondence principle, Kepler's orbits, Eq. (13), must be recovered 
in the limit e — > 0, so that 3eA = e is the eccentricity of Newtonian mechanics. To first 
order in e, Eqs. (18)-(20) become 

f c wr c (l-3e), (21) 
gpae(l + 3e), (22) 
« ^ 1 - 3e, (23) 
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and the relativistic orbits in this limit may be expressed concisely as 

rc(1 - 3e) « 1 + e (l + 3e) cos (1 - 3e)y>. (24) 
r 

Equations (21) and (22) are the relativistic radius of circular orbit and relativistic eccentric- 
ity, respectively, and are discussed in more detail in Sec. III. 

The approximate orbit equation in Eq. (24) predicts a shift in the perihelion through an 
angle 

Aip = 2n(R- 1 - 1) « 2vr(3e) (25) 

per revolution. This first-order prediction is in agreement with well known results [1, 2, 3, 

4, 5, 6, 7, 8, 9, 10, 15, 23, 24, 25, 26, 27, 28, 29], and agrees with the observed precession 

of perihelia of the inner planets [1, 2, 3, 4, 5, 6, 9, 15, 23, 30, 31, 32, 33]. (Precession 

due to general relativity is illustrated in Fig. 1.) As with other first-order calculations, 

this result may be compared to observations, assuming that the relativistic and Keplerian 

angular momenta are approximately equal [1, 2, 3, 4, 5, 6, 7, 8, 9, 15]. For a Keplerian orbit 

[20, 21, 22] i 2 = GMa(l - e 2 ), where G = 6.670 x lO" 11 m 3 kg" 1 s~ 2 , M = 1.989 x 10 30 kg 

is the mass of the Sun, and a and e are the semi-major axis and eccentricity of the orbit, 

respectively. Therefore, the relativistic correction defined after Eq. (11), 

3GM 
c z a[l — e z ) 

is largest for planets closest to the Sun and for planets with very eccentric orbits. For 
Mercury [34] a = 5.79 x 10 10 m and e = 0.2056, so that 3e « 7.97 x 10~ 8 . (The speed of 
light is taken to be c 2 = 8.987554 x 10 16 m 2 /s 2 .) According to Eq. (25), Mercury precesses 
through an angle 

A<^~ — = 5.02 x 10' 7 rad 27 

c 2 a(l — e 2 ) 

per revolution. This angle is very small and is usually expressed cumulatively in arc seconds 

per century. The orbital period of Mercury is 0.24085 terrestrial years, so that 

100 yr 360 x 60 x 60 A 

A$ = — x x Atp (28a) 

0.24085 yr 2vr r K ' 

?s 43.0 arcsec/century. (28b) 

Historically, this contribution to the precession of perihelion of Mercury's orbit precisely 
accounted for the observed discrepancy, serving as the first triumph of the general theory of 
relativity [11, 12, 13, 14]. 



III. DISCUSSION 



The approximate orbit equation in Eq. (24) provides small corrections to Kepler's orbits 
due to general relativity. [Compare Eqs. (24) and (13).] A systematic verification may be 
carried out by substituting Eq. (24) into Eq. (11), keeping terms of orders e, e, and ee only. 
The justification for discarding the term nonlinear in the eccentricity is the correspondence 
principle. Arguments concerning which terms to discard based only on direct comparisons 
of relative magnitudes of higher-order and lower-order terms lead to contradictions (as may 
be verified by the reader). Instead, the domain of validity is expressed by subjecting the 
solution Eq. (24) to the condition 

- - 1 « 1 (29) 
r 

for the smallest value of r. Evaluating the equation of the orbit Eq. (24) at the perihelion 
(r p ) results in 

r c _ l + e(l + 3e) 

r P ~ 1-36 • (30) 
The substitution of this result into Eq. (29) results in the domain of validity: 

e(l + 3e) +2(3e) < 1. (31) 

Therefore, the relativistic eccentricity e = e(l + 3e) <C 1, and Eq. (24) is limited to describing 
relativistic corrections to near-circular (Keplerian) orbits. Also, the relativistic correction 
2(3e) <C 1, and thus the equation of orbit Eq. (24) is valid only for small relativistic correc- 
tions. 

Equation (24) is also consistent with a reduced radius of the circular orbit, as derived 
from the Schwarzschild effective potential: 

V ^~— + ^-^- (32) 

This effective potential is customarily defined [3, 5, 7, 8, 9, 10, 15, 35] after making the obser- 
vation that Eq. (7) reduces to the Newtonian result as r s — > (c — > oo). [The Schwarzschild 
effective potential in Eq. (32) is compared to that derived from Newtonian mechanics in 
Fig. 2.] If we minimize V e s with respect to r, we obtain the radius of the circular orbit, 

Rc = \r c + \r c ^\ - We « r c (l - 3e), (33) 
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in agreement with Eq. (21) to first order in e. An additional characteristic of relativistic 
orbits is that of increased eccentricity Eq. (22). Equation (24) predicts that relativistic 
orbits will have increased eccentricity compared to Keplerian orbits. The correction to the 
eccentricity is of the same order as that for the relativistic contribution to perihelic precession 
and the reduced radius of the circular orbit. 

IV. CONCLUSION 

The relativistic central-mass problem has been cast into finding a solution that is not 
very different from a Newtonian circular orbit, the result of which is a Keplerian orbit with 
small relativistic corrections. The resulting first-order relativistic contribution to the peri- 
helic precession in Eq. (25) agrees with established calculations. This effect is one of several 
relativistic effects which arise as first-order corrections to the familiar orbit equation describ- 
ing Keplerian orbits in Eq. (13). The approximate orbit equation is limited to describing 
small relativistic corrections to near-circular Keplerian orbits, as expressed by the domain of 
validity in Eq. (31). Comparisons with observations are made with the standard assumption 
that the relativistic and Keplerian angular momenta are approximately equal. 

Compared to the usual perturbative approach, this derivation demonstrates several de- 
liberate steps toward a Keplerian limit prior to solving the nonlinear differential equation 
(10) arising from the Schwarzschild line element. The linearization procedure preceding 
Eq. (14) and appeal to the correspondence principle [to identify the Newtonian eccentricity 
in Eq. (19)] makes for a very physical approach toward the boundary between relativistic 
and Keplerian orbits. The resulting approximate orbit equation in Eq. (24) lends itself to 
easy comparison with the familiar Keplerian orbits in Eq. (13) and displays several rela- 
tivistic effects. This simple solution arises as a result of pursuing a narrower problem from 
the beginning, and explicitly expressing the domain of validity in Eq. (31). This derivation 
is approachable to undergraduate students, and the resulting orbit equation provides quali- 
tative and quantitative understanding of the corrections to Keplerian orbits due to general 
relativity. 
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FIG. 1: Relativistic orbit in a Keplerian limit (solid line) as described by Eq. (24), compared to 
a corresponding Keplerian orbit (dashed line), Eq. (13). The precession of the perihelion is one 
characteristic due to general relativity, and is illustrated here for < <p < 20ir. The eccentricity is 
chosen to be e = 0.2 for both the relativistic and Keplerian orbits. This characteristic is exaggerated 
by the choice of relativistic correction parameter (e = 0.1) for purposes of illustration. Precession 
is present for smaller (non-zero) reasonably chosen values of e as well. 
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FIG. 2: The effective potential in Eq. (32) derived from the Schwarzschild geometry (solid line) 
compared to that derived from Newtonian mechanics (dashed line). The vertical dotted lines 
identify the radii of circular orbits, R c and r c , as calculated using the Schwarzschild geometry and 
Newtonian mechanics, respectively. The Schwarzschild geometry predicts a smaller radius of the 
circular orbit than that predicted by Newtonian mechanics [Eq. (33)]. The value e = 0.08 is chosen 
for purposes of illustration. 
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